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Abstract 

We have studied the interference of a variable number of independently created rriF = micro- 
condensates in a C02-laser optical lattice. The observed average interference contrast decreases 
with condensate number N. Our experimental results agree well with the predictions of a random 
walk model. While the exact result can be given in terms of Kluyver's formula, for a large number 
of sources a 1/yN scaling of the average fringe contrast is obtained. This scaling law is found to 
be of more general applicability when quantifying the decay of coherence of an ensemble with iV 
independently phased sources. 
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Interference experiments are valuable tools for the study of coherence of atomic sources Q, 
While for an incoherent ensemble the number of independently phased objects equals the 
atom number, in a Bose-Einstein condensate atoms macroscopically populate a one-particle 
quantum state and the whole condensate is a single coherent atom source. Of large interest 
is the intermediate region of partial coherence. We shall here be interested in the degree of 
coherence of an ensemble of N independent quantum objects. For an array of periodically 
spaced atom sources, each of which be coherent, Ashhab predicted a l/y/N scaling of the 
average fringe visibility with the number of sources N in a far field interference experiment 
0| . For other theoretical works, see U Experimentally, Hadzibabic et al. had observed 
the interference from 30 independent condensates and obtained the surprisingly high average 
fringe contrast of 34%, which was indeed comparable to the results observed in the early 
two-condensates interference experiment of [lj]. The experimental result of Hadzibabic et al. 
was in agreement with a Monte Carlo simulation. 

We here experimentally investigate the interference of a variable number of randomly 
phased atomic sources. To verify the predicted 1/yN scaling of the fringe visibility as 
a function of the number N of sources, an array of magnetic field-insensitive (mj? = 0) 
Bose-Einstein condensates is created by all optical cooling and trapping techniques. In the 
sites of a one-dimensional mesoscopic optical lattice potential, rubidium atoms are indepen- 
dently cooled evaporatively to quantum degeneracy. After releasing the microcondensates, 
an interference pattern is created. The observed average fringe visibility decreases with the 
number of coherent sources over which the atoms are distributed. Our results well follow 
the predictions of an exact theoretical model based on a random walk in the complex plane. 
In the limit of large number of condensates N, the expected 1 / \/N scaling of the average 
fringe visibility characteristic for the central limit result of the random walk is experimentally 
reproduced. 

The observed variation of fringe visibility with the number of sources can intuitively be 
understood as follows. When two independent condensates with equal number of atoms are 
overlapped, upon measurement the condensates are projected onto a state with well-defined 
relative phase. This phase is random in each realization of the experiment, but in all cases 
a density modulation with 100 % contrast occurs ODD . For three interfering condensates 
now two relative phases come into play, and in general the phase difference between any two 
of the condensates differs from that between either one of them and the third condensate. 
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One still expects a periodic density modulation, but with reduced contrast. For a larger 
number of independent sources, the degree of randomness increases. The fringe contrast for 
an array of N coherent sources can be calculated with a random-walk model in which the 
number of steps equals the number independent phases, i.e. N — 1 jjj. 

We wish to point out that the predicted 1/vN scaling of the fringe visibility with the 
number of sources N in the lattice is actually a far more general result when considering 
the degree of coherence of atomic ensembles. As a textbook-like example, let us consider 
a sample with N independent monochomatic sources (single atoms in a classical model or 
independent condensates) with total wavefunction $(x) oc Yln=i e^ k x+6,?l - ) , where the phases 
6 n (n = 1, . . . , N) are randomly distributed. To test for the coherence of this ensemble, 
we shall beat the sample with a coherent local oscillator wave (e.g. an atom laser) with 
wavefunction <&i (x) oc e iklo ' x . This yields a spatial interference signal |$+$i | 2 with intensity 
modulation: $* $ + c.c. = 2 Re (e i ( k_kl °)' x ^ n =i e%6n "\ After evaluating the average value of 
this random walk problem with N steps, one finds that the modulation of the fringe signal 
scales as V~N, which compares to the value of N that is obtained with a coherent atomic 
ensemble. The visibility of the interference signal, a quantity commonly used when quantify 
coherence in the optics literature Q], decays with a l/y/N scaling law, i.e. relatively slowly 
with the number of sources. 

In our experimental approach, atoms are evaporated to quantum degeneracy in the sites 
of a one- dimensional periodic potential created by an optical standing wave generated by 
C02-laser whose wavelength is near 10.6 fim. This laser frequency is roughly a factor 14 be- 
low the lowest electronic resonance frequency of the rubidium atoms used in our experiment. 
Due to the static atomic polarizability, atoms are pulled into the antinodes of the standing 
wave. Adjacent sites are spaced by d — Aco 2 /2 — 5.3 fim, which is sufficiently far apart that 
tunneling between sites is completely negligible. A variable number of sites is initially occu- 
pied with thermal atoms, which are subsequently cooled to quantum degeneracy by direct 
evaporation in the optical trap. Our experiment benefits of recent developments in the "all- 
optical" generation of Bose-Einstein condensates in optical dipole potentials |lfl, 11 , 12, 3], 



which we here extend to a lattice geometry. Besides offering the flexibility to directly evap- 
orate into the periodic atom potential, in such traps magnetic-field insensitive condensates 
can be produced Q , which is a clear benefit in studies of coherence properties. 

Let us theoretically discuss the expected interference signal after releasing the array 
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of N (with N variable) independent microcondensates in the far field. This treatment is 
inspired by earlier calculations considering the interference of independent optical 14] and 
atomic sources. In our model, we consider the N independent atomic Bose-Einstein 

condensates with same trap parameters in all sites before the free expansion. The total 
wavefunction at this time can be written as 

N 

$(x) =Y,a n e ie "e-( x - x ^^ (1) 

n=l 

where a n e l9n denotes the amplitude of the n-th condensate wavefunction, which is assumed 
to be a gaussian function with ground state size I = \Jh/ (muj). The microcondensates 
are centered at the C02-laser lattice antinodes, corresponding to the positions x n = n ■ d, 
where n — 1, . . . , N. Because the microcondensates are independent, the phase factors 9 n are 
randomly distributed . We can immediately derive the wave function in a momentum 

space picture by performing a Fourier transform. When expanding the atomic clouds 
and recording the interference pattern with a time of flight measurement in the far field, we 
actually perform a measurement of |$(p)| 2 , which can be written as 



2 r N-l 

2 I 



A n cos(ndp/ h + <f n ) 



n=l 



e~^' h \ (2) 



In this equation, Aq = Na 2 , tpo — 0. For n other than zero the (real) amplitudes A n and 
the corresponding phase angles (p n are given by jjj A n e l(fn = 2a 2 Y2^=i 1 ^ 0q ~ eq+n \ if we for 
the sake of simplicity assume equal atom population for all sites (i.e. a n = a for all n). 

In our experiment, we only spatially resolve the fringe pattern with largest spatial mod- 
ulation, i.e. the interference between nearest neighbours. The expected visibility of this 
modulation is given by V = (J max — -^min)/(-^max + -^min) = Ai/A , which we can write more 
simply as V = 2S/N, where S = \z\ = \J2q=i z q\ = \J2q=i e lA9q \ and (pi = ip equal 
the phase angle of this sum. The phase angles A8 q are randomly distributed in the range 
< A6 q < 2tt. For each realization of the experiment, the fringe contrast is determined by 
the modulus of a sum of unit vectors z q which have random directions in the complex plane. 
Further, the phase angle of the total vector gives the position of the fringe pattern. In our 
model, the problem of calculating the fringe contrast and position of an array of coherent, 
independent sources thus can be transferred to the solving of a random walk problem in the 
complex plane. This is a typical two-dimensional Pearson's random walk Q]. Fig. la 
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gives a graphical representation of the situation. In general, the expected probability density 



to arrive at a distance S from the origin is given by Kluyver's formula 



12|: 



p( s > n ~ !) = A rrr J o [Hp)] N ' 1 HS P )pdp (3) 

where Jo is the 0-th Bessel function. For large N, the Bessel functions can be asymptoti- 
cally expanded, and a Gaussian probability distribution is obtained, corresponding to the 
celebrated central limit result |lfij |. In this approximation, the expected average visibility 
approaches the previously derived value of y/n/(N — 1) ~ ^ir/N, when N 3> 1 0. Note 
that N — 1 equals the number of relative phases between adjacent sites and therefore the 
number of steps of the random walk. In the general case, and especially for a small number 
of condensates, one has to fully solve Eq. |3] which is difficult due to the oscillatory character 
of the Bessel functions. We have numerically evaluated the probability density by expand- 
ing Kluyver's formula in a Bessel- Fourier series Q|. As an independent verification of our 
calculation, we performed a Monte Carlo analysis of the average visibility assigning random 
values to the arguments A6 q . Fig. lb shows the theoretical expected fringe average visibility 
as a function of condensate number derived using the exact calculation (dots), the central 
limit result (dashed line) and the simple ^ir/N scaling function (dashed-dotted line). 

An array of independent Bose-Einstein condensates with variable number of sites is pro- 
duced using a modification of our previous setup [3, Q . In an ultrahigh vacuum chamber, 
mid-infrared radiation near 10.6 pm for atom trapping is focused to a beam waist (1/e 2 
radius) in the range 25-40 pm. With a second lens and a retroreflection mirror, a standing 
wave is generated. The size of the beam waist in the trapping region could be varied with a 
telescope. This gave us control over the number of populated trapping sites. For a loading 
of atoms, we switched to a running wave geometry by slightly misadjusting the CCvlaser 
beam backreflection mirror away from a perfect retroreflection. This mirror was mounted 
on a piezo, which allowed for an electronic control of the degree of alignment during the 
course of the experiment. Cold rubidium ( 87 Rb) atoms from a magneto-optical trap (MOT) 
were loaded into the running wave geometry. After this transfer, 10 6 atoms, populating 
the lower hyperfine ground state (F = 1, mp = 0, ±1) are left in the optical trap at a 
temperature near 100 pK. To cool the trapped atomic cloud, the C02-laser beam power is 
acousto-optically ramped down to induce forced evaporative cooling. The highest energetic 
atoms leave the purely optical trap, and the remaining ones thermalize to lower tempera- 
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tures. The total evaporation stage lasts about 10 seconds, during which the mid-infrared 
beam power is smoothly reduced from 30 W to a typical final value of 40-50 mW. Through- 
out the evaporation stage, the MOT quadrupole field with 10 G/cm field gradient is left on. 
This gradient is sufficiently strong to remove atoms in field-sensitive spin projections in this 
running wave geometry, as the trapping force along the beam axis is here relatively weak [19] . 
While the initial phase of this evaporation is still performed in the running wave geometry 
(this limits the number of sites over which atoms are distributed), during the course of the 
evaporative cooling we slowly switch to a standing wave (lattice) geometry by servoing the 
piezo-mounted mirror correspondingly. The ID lattice geometry is fully aligned at a time 
1-2 seconds before the onset of quantum degeneracy (corresponding to atomic temperatures 
of a factor 2 above the transition temperature of about T c ~ 200 nK). This ensures that the 
microcondensates are formed independently. 

At the end of the evaporation state, an array of mp = microcondensates is created. 
The total number of atoms in the optical lattice is about 7000, and the number of condensed 
lattice sites can be varied from typically 5 to 35 by choosing different beam waists of the 
lattice beams and - in a smaller range - also allowing for tiny residual misalignments of the 
lattice beams at given beam diameter. For all our measurements, the estimated tunneling 
time between sites is above 10 18 s. This assures that the microcondensates are truly inde- 
pendent. The residual sensitivity of the mp = condensates to stray magnetic fields due 
to the second order Doppler shift is near 14fK/mG 2 . For a typical extension of the lattice 
of 100 fim, the variation of the condensate phases due to magnetic field inhomogenities is 
0.15 Hz at an estimated field gradient of 50mG/cm. On the other hand, we expect that 
differential mean field shifts due to variation in the atom number do cause a non-negligible 
differential variation of the condensate phase with time. 

After creation of a variable number of microcondensates in the lattice, the C02-laser 
radiation is extinguished, and the atomic clouds are allowed to freely expand and fall in 
the earth's gravitational field. The condensed atomic clouds originating from adjacent sites 
overlap as soon as the free expansion time t exceeds 0.5-1.5 ms, depending on the interaction 
energy for the particular atom population per site, which in turn depends on N. We then 
expect to observe an interference pattern with spatial period of A t h — ht/ (md), despite each 
individual atom source having random phase. Figure 2a shows a typical observed interference 
pattern for N = 20 interfering coherent atomic clouds. The free expansion time here was 
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t = 15 ms, after which the shown absorption image was recorded. A horizontal profile of this 
image is shown in Fig. 2b, where the dashed line with dots shows the experimental data. 
For an analysis of our data, the fringe patterns (i.e. the absorption image profiles) were 
fitted with the function 



where A exp and (p exp denote the fitted fringe spacing and the phase angle of the interference 
pattern. Gaussian envelopes with independent widths and positions are here included for 
the background and fringe signals respectively. The fringe visibility V exp of the pattern 
of Fig. 2a, as derived from the fit, was 31.5 %. The experimental fringe spacing here is 
A exp — 13.0 ± 0.3 /mi, which agrees well with the expected spacing A t h — 12.97/xm. We 
have studied the variation of fringe contrast and phase angle for different realizations of the 
experiment. Fig. 2c shows these parameters in a polar plot, where each dot corresponds to 
the result of an individual realization. Phase angles and fringe contrasts here vary randomly 
from measurement to measurement. Moreover, no preference for a certain value of phase 
angle is visible. Note the close analogy of this polar image to Fig. la. The vector z = Se lip as 
a complex signal modulation (where S = y V) represents the endpoint of a random walk in 
the complex plane, and can be directly experimentally determined from fringe contrast and 
phase of a far field interference pattern. Series of fringe patterns were recorded for different 
numbers N of populated lattice sites. 

Fig. 3 shows the variation of the average fringe visibility for a series of measurement as a 
function of N. The average number of populated sites here was determined by recording the 
length of the lattice, as measured with a time of flight image at zero expansion time. Our 
experimental data clearly shows a decrease of the fringe contrast with condensate number, 
which agrees with the exact random walk solution valid for arbitrary iV (see Eq. |3J) and 
within our experimental uncertainties also with the central limit result. For large number of 
independent sources, the fringe contrast averages out. However, this averaging effect occurs 
relatively slowly, i.e. it takes a large number of sources to let the fringe pattern resemble 
that of an incoherent atomic sample. The data shown in Fig. 3 has been fitted with the 
(interpolated) result of the exact calculation, as shown by the solid line. To account for our 
experimental imaging resolution, we have multiplied the theoretical result with a constant 




(4) 
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value, which was the only free parameter in the fit. The resolution of 6.3 /zm obtained in 
this way agrees well with an independent, direct measure of the optical resolution of 6.0 /im. 
The dashed line gives the corresponding visibilities using the simple ^n/(N — 1) result of 
the central limit approximation, which was multiplied with the same constant factor. 

To conclude, we have studied the interference of an array of independently generated 
mp = microcondensates in an optical lattice for a variable number N of populated sites. 
The observed average interference visibility decreases with source number due the increased 
number of (upon measurement) randomly distributed relative phases, and well follows the 
predicted 1/yN scaling. 

For the future, we anticipate that the observed relatively slow decay of coherence with 
source number will allow one to observe related interference effects with small samples of 
thermal atoms. A further direction for future research would be to experimentally verify the 
intrinsic (i.e. quantum-mechanical) randomness of the relative phases of the atom sources. 
Due to the use of mp = condensates, differential phase shifts due to stray magnetic fields 
are estimated to contribute clearly below a radian. Our current experiment, differential 
mean field phase shifts due to the in general slightly different atom numbers per site during 
the approximately 1-2 seconds time period between the onset of quantum degeneracy and 
the release of the microcondensates are however estimated to still accumulate statistical 
variations of the phase p of several tens of radians from shot to shot. If e.g. measurements 
with constant atom number per site could be postselected, which clearly requires small 
absolute atom numbers, the intrinsic randomness of the condensates relative phase might be 
experimentally tested. Finally, we anticipate that our experimental scheme holds promise 
for novel types of lattice atom lasers. 
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FIG. 1: (a) Random walk of coefficients z q = e lAdq , with the phase angles < A9 q < 2ir arbitrary 
in the complex plane. The vector z is the sum of N — 1 independent coefficients z q (I < q < N — 1) 
and determines both the fringe visibility V = 2S/N with S = \z\ and its phase angle tp. (b) 
Predicted average fringe visibility for the far field interference pattern of N coherent atomic sources 
(dots). Result of the central limit approximation: V = yJn/(N — 1) (dashed line) and the simple 
V = \/tt/N scaling function (dashed-dotted line). 

FIG. 2: (a) False-color absorption image of an interference pattern of 20 independently generated 
mj? = microcondensates released from a C02-laser optical lattice. The image was recorded after 
a free expansion time of 15 ms. The field of view is 350 /mi x 140 /im. (b) Horizontal profile of 
the image averaged over a vertical region of 35 /im: Experimental data (solid line with dots) and 
fitted fringe pattern (dashed line), (c) The fringe visibilities and phase angles of fringe patterns 
arising from the interference of 20 independent Bose-Einstein condensates, as derived from fits as in 
Fig. 2b, have been drawn (dots) in a polar diagram for 130 different realizations of the experiment. 

FIG. 3: Average fringe visibilities of far field interference patterns as a function of number of 
interfering coherent atomic sources. The experimental data sets were recorded with beam waists: 
24.3 /im (squares), 30.0 /im (dots), and 40.2 /im (circles). The data has been fitted with the (inter- 
polated) theoretical fringe contrast multiplied by a constant factor, which is left as a free parameter 
to account for the finite imaging resolution (solid line). The dashed line gives the corresponding 
result in the central limit approximation when assuming the same imaging resolution. 
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